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ABSTRACT 

Numerical studies of ion-dipole collisions are extended to calculating 
the effects of vibrational degrees of freedom on capture cross sections and 
formation of ion-molecule complexes. Two different dipole moment varia- 
tions about the equilibrium intemuclear separation of the oscillator are 
studied. The capture cross section is sensitive to the oscillator model. 
Computer-made motion pictures reveal a change in mode of the oscillation 
due to the ion -dipole interaction. The oscillator becomes constrained so 
it cannot relax to its equilibrium separation. This shift in oscillator 
potential barrier is denoted as "classical tunneling" of the oscillator. 
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INTRODUCTION 

Capture cross sections and collision complex lifetimes have been cal- 
culated numerically for ion-polar molecule collisions." 1 ' ^ The energy 
transfer in these collisions has been studied using computer -plotter tech- 
niques.^ ^ (The polar target molecule is represented by a rotating dipole 
imbedded in a polarizable sphere. The ion is made to reflect at an inter- 
molecular distance r c to simulate the hard repulsive cores of the ion 
and molecule.) Capture cross sections have been calculated to set an up- 
per limit to ion-molecule reaction cross sections and to compare the pre- 
dicted ion energy dependence with mass spectrometric results.^ ® The 

3 

permanent dipole introduces multiple reflection behavior corresponding to 
the formation of long-lived ion-molecule collision complexes. These com- 
plexes may be of importance as clustering sites in gas discharges or ir- 
radiated gases. (Such "captures" have calculated collision lifetimes as 
large as several hundred times a single reflection period. ) Motion pic- 
tures have also been made with the computer to study hindered rotation of 
the dipole by the incident ion.^ Hindering may favor a specific chemical 
reaction. 

Lifetimes for multiple-reflection collisions have previously been 
calculated including rotational and translational degrees of freedom for 
the ion-molecule collision pair. This paper reports on numerical studies 
of vibrational effects on the capture cross section, the probability of 
multiple reflections, and collision lifetimes. The vibrational effects 
are included by approximating the polar target as a (rotating) harmonic 
oscillator whose dipole moment is a function of internuclear separation 
along the molecular axis. The oscillator motion is followed in a 
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coordinate system whose origin is fixed, at the center-of-mass (CM) of the 

target molecule. The CM of the oscillator is at the center-of-mass of the 

dipole and the oscillator separation is much less than the ion-molecule 

separation. (This is consistent with the ideal dipole approximation.) 

COLLISION MODEL AND INITIAL CONDITIONS 

At high frequencies (i.e., on the order of molecular vibrations and 

fast rotations) the dipole moment of a polar molecule is not constant but 

a, function of its internuclear separation. The separation of interest 

9 

is usually that along a particular interatomic bond in the molecule. The 

variation in dipole moment y is assumed gaussian for these studies (i.e., 
_ 2 2 

y = y^e ), where y^ is the "static" dipole moment, a is the damping 

constant, and x is the reduced bond distance (r - r where r is the 

e 

instantaneous separation and r g is the equilibrium bond separation). 

Such a variation is a reasonably representative approximation for diatomic 
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and polyatomic molecules of interest. To test the sensitivity of the re- 
sults to the shape of the y(x) variation, the damping coefficient a was 

o 

selected so that y = 2y /3 at x = ±0.05 A for oscillator set (A). The 

o 

maximum x value for the free oscillator assigned was 0.05 A (0.005 nm) 

o 

in all cases. For oscillator set (B) y = 0.02 • y Q at x = ±0.05 A. The 

latter approximation is probably more realistic for polar targets. The 

shapes of the two gaussian variations are shown in figure 1. The period 

of change in y is Wn. /k , where k and m are the oscillator 

o r o o o o 

force constant and reduced mass. The values of k and m for the os- 

o o 

dilator set the maximum values of x and x at infinite separation for 

fixed E . The m value is not the same as the mass of the target mole- 
o o 

cule. Its value was chosen so as to reduce computer time. The energy 
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values E q chosen were 0.3 6 eV and 0.03 6 eV corresponding to stiff and 

soft (nearly thermal) oscillations. However, only 50 cases were studied 

for the former because computer time was prohibitive. A k Q value of 

4.8x10^ dynes cm ^ (0.48 Newtons/cm) was assigned for these calculations 

with an m Q value equal to ten times the mass of the hydrogen atom. 

This value is one-tenth the ground state value for HC1. Integration 

for trajectories of all oscillators causes extended computer running 

times (minimum of 5 min. for v = 5x10^ cm .sec ^). This is because the 

-2 

oscillation period is only 10 of the characteristic rotational period 

where = 2 ^/l/2E^; E R is the rotational energy and I is the 
moment of inertia. The oscillator motion somewhat separates out from the 
rotational and translational motion; however, it will be shown that the 
nature of the variation is important in determining the ion turning 
points. Individual collisions with stiff targets require as much as 10 
minutes to one-half hour computer time which is clearly prohibitive if 
hundreds of cases must be studied. 

Most of the collisions studied involved the soft oscillators with 
E q = 0.036 eV. The total energy of the free oscillator is 

E = jj- m x 2 + ^-k x 2 (l) 

o 2 0 2 o 

where all quantities have been previously defined. The maximum displace- 
ment of the free oscillator at large distances from the ion is given by 
(x n>o = hfih i.e., the separation where the oscillator energy is 
purely potential with no interaction. The initial separations x were 

chosen randomly in the interval |x| < ( x m ) Q * Approximately 400 collisions 

o 

have been studied for impact parameters B from 3 to 12 A. These b 
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values correspond roughly to critical impact parameters for capture with 

1-3 

CO and CH^CN targets. The oscillator energy is constrained to remain 
positive throughout interaction but is allowed to decrease in value. 

Such behavior is, of course, not possible for a quantum oscillator in 
the ground state which must have energy 7 j-hv Q . Thus the observed 
"cooling" of such oscillators applies only to molecules which are ini- 
tially in excited vibrational states. It should also be noted that the 
classical oscillator approximation is least adequate in describing vari- 
ations in energy of stiff oscillators since quantum gaps are relatively 
large. The differential equation of motion for the oscillator is 

X = - 2a^(ye cos y/R^ + k Q )x/m o (2) 

The first term is an interaction term whereas the second term is a con- 
ventional restoring force term for the oscillator. The relative magni- 
tude of these terms is important in determining dipole behavior. 

RESULTS 

The probability of multiple reflection f D for oscillating polar 
targets is composed with results for purely rotating polar targets f 

n 

2 

in table I. The results of set A are for an a value which makes 

o 

y = 2y Q /3 at x = ±.05 A. Results are presented for both CO and CH^CN. 

2 3 °-2 

Results for set (B) are for an a value (l.6xl0 J A ) which makes 

o 

y ~ 0.02y o at x = ±.05 A and include only CH^CN results. The effect 
of oscillations clearly depends on impact parameter, however, it is dif- 
ficult to detect general trends in the results. The f value is uni- 

K 

formly higher than f_, for set (b) at several b values but drops to 

O 

zero at b = 11 A. The variation y(x) for the oscillator is important; 
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this can be seen from the difference between identical random number sets 
2 

for different a values at roughly the same b values. 

Further, there is a dramatic effect on the capture cross section 

2 

as demonstrated by the plots of CL vs. b for CH_CN in Fig. 2. There 

n d 

2 

are fewer points for the vibrating targets (a for y ~ 0.02y Q at 
o 

x = ±0.05 A and only 25 cases per point (compared to 50 cases) because of 

the prohibitive computer running times. The results agree at low impact 

o 

parameters b =6 and 8 A. However, there are marked differences between 

o o 

results for b = 10 A as well as for b = 11 A. The differences between 
values are well outside the rough 10# variance in cross section be- 
tween 25 and 50 case sets of trajectories. The oscillator studies were 

o 

not done for b > 11 A since the C R value approaches zero. The low- 
ered C_ values are for "heated" oscillator collisions so the result is 

Ja 

meaningful for real molecules in the limits of the classical interaction 

approximation. The effect of the oscillator on capture cross section will 

be much less for CO than for CH^CN since the polarizability term determines 

4-5 

the capture cross section for CO. 

The collision lifetimes are relatively insensitive to vibrations for 
a given species and the average for multiple reflection collisions is 

~2 to 3x10 ^ sec for CH^CH targets. It should be noted that the trajec- 
tories were cut off at a maximum of 10 reflections so that no meaningful 
maximum reflection time could be calculated. Results of previous studies 
indicate that the average number of (multiple) reflections is 15 to 25 for 
certain b values.^ 

Certain features of the energy exchange between the translational, 
rotational, and vibrational degrees of freedom are shown in the plots of 
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Fig. 3. The modulation of the interaction due to the variation of y 
with x is almost undetectable because the vibrational time scale is so 
small. 

The modulation of the oscillator separation x via the ion-dipole 

o 

interaction is evident in the region R ~ 5 A as shown in fig. 4(a) for 
a CH^CN target. The corresponding variation in oscillator energy is 
shown in fig. 4(b). 

Multiple reflection behavior is shown in fig. 5(a) and. for a repre- 
sentation collision. The hindered rotation accompanying multiple reflec- 
tion is shown in fig. 5(b). The symmetrical nature of the net spiraling 
(due to multiple reflections) is somewhat distorted by the oscillator. 
However, the general features of this phenomenon are 'unchanged by the 

oscillator as is hindering of the dipole. The hindering of the CH^CN tar- 

o 

get occurs at large ion-molecule separations (»7 -8 A) as can be seen from 
the decreasing envelope of the T plot in fig. 5(b). The characteristic ion 

O 

turning points of the ion-molecule complex are large, ~ 10 to 20 A for 

ch 3 cn. 

COMPUTER MOTION PICTURES - "CLASSICAL TUNNELING" 

Vibrational motion has been added to the ion-dipole movies which 
were described in refs. 5 and 10. The modulation of the oscillator which 
is suggested in the time history plots is confirmed in the movies . This 
modulation is of course dependent on the value of cos y, i.e. , the ion- 
dipole orientation angle which depends on their relative motion. 

In viewing the movies it was noted that there is a change in the mode 
of the oscillation in several instances. This change in mode occurs when 
the oscillator acceleration goes through zero; that is, when the restoring 
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force is balanced by the interaction term. This constraint of the oscilla- 
tor so that it cannot relax to its equilibrium is somewhat analogous to 
quantum mechanical tunneling through a potential barrier so we have chosen 
to call it "classical tunneling". The time history plots of oscillator 
coordinate x and velocity x in figure 6 show that X stays near zero 
for a significant fraction of an oscillation period (where x is nearly 
constant) . 

CONCLUDING REMARKS 

Preliminary studies of vibrational effects on multiple reflections, 
collision lifetimes, and capture cross sections in ion-dipole collisions 
have been completed. For thermal oscillators (E q = 0.03 6 eV) the effect 
of the vibrator depends on impact parameter as well as target species . 

The fraction of multiple-reflection cases is a complicated function of 
the oscillator model, particularly the shape of the dipole moment varia- 
tion p(x). The oscillator can have a dramatic effect on the capture 
ratio C^; this lowering of the capture cross section is particularly 
sensitive to the damping constant a. The inclusion of additional degrees 
of freedom, particularly if coupled, might markedly increase the collision 
complex lifetimes. Vibrator studies must be extended to a variation in 
parameters before more general conclusions can be drawn about the role of 
the oscillator. "Classical tunnelling" of the oscillator via the ion- 
dipole interaction suggests a mechanism for rearrangement and/or dissocia- 
tion which operates without the force constant changing. Quantitative 
conclusions regarding this and related phenomena require further study. 
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TABLE I - MULTIPLE REFLECTION RESULTS FOR HARMONIC OSCILLATOR STUDY 


FOR ION VELOCITY OF 5x10^ CM SEC 1 , 25 COLLISIONS 


STUDIED PER IMPACT PARAMETER 


Polar 
molecule 
set (A) 

Impact 

parameter, 

0 

A 

( IQ-10 nm ) 

Fraction of 
collisions with 
mult. ref. 

f R 

(without oscillator) 

Fraction of 

collisions with 

mult. ref. 

f ' 

R 

(with oscillator) 

CO 

3.0 

0.24 

0.16 

CO 


.20 

.16 

CO 

5- 

.2k 

.16 

ch 3 cn 

5. 

ro 

.20 

CH 3 CN 

6. 

.16 

.08 

CH 3 CN 

T. 

.36 

.48 

CH 3 CN 

9- 1 

.28 

.36 

CH CN 

11. 

.20 

.20 

set (B) 




CH 3 CN 

6.0 

0.16 

0.60 

CH 3 CN 

8. 

.24 

.12 

CH 3 CN 

10. 

.24 

.32 

CH 0 CN 

11. 

.20 

.0 

3 j 
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(A) a 2 - 160 A~ 2 . (B) a 2 = 1. 6X10 3 A' 2 . 


Figure 1. - Dipole moment variation for polar target oscillators studied (sets A and B), |x m | 0 = 0.05 A ; 
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Figure 2. - Comparison of capture ratio variation with impact parameter for oscillating and pure rotating polar 
CH 3 CN targets from set (B). Results are for an initial ion velocity of 5xl0 4 cm sec " 1 where the cross section 

o r =jr /' oc Cn(b 2 )d(b 2 ). 



ROTATIONAL ENERGY SCALED VELOCITY SCALED TO INITIAL 

TO KTn = 0. 026 eV VALUE OF 5. 5X10 4 CM SEC' 1 




Ion -morl ecu !e separation, A 


Figure 3. - Variation of ion velocity and polar molecule rotational energy 
during methyl cyanide-parent ion multiple reflection capture collision. 
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(A) PROJECTIONS OF ION TRACE 
IN CARTESIAN XY PLANE OF 
CENTER -OF-MASS SYSTEM. 



(B> ION-DIPOLE ORIENTATION ANGLE. 


Figure 5. - Evidence for multiple reflection and hindered rotation behavior 
in CH 3 CN ion-dipole collisions with oscillating target. 



Figure 6. - Variations of oscillator coordinates and velocity for ion-dipole collision demon- 
strating "classical tunneling". 
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